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Abstract The polarimetric and photometric variability of Wolf-Rayet (WR) stars as caused 
by clumps in the winds, is revisited. In the model which is improved from Li et al. 2000, 
the radial expansion of the thickness is accounted for, but we retain the dependence on the (3 
O ' velocity law, stellar occultation effects. We again search for parameters that can yield results 

consistent with observations in regards to the mean polarization p, the ratio TZ = tTp/iTphot 
^ ■ of polarimetric to photometric variability, and the volume filling factor /y. Clump generation 

and spatial distribution are randomized by the Monte Carlo method so as to produce clumps 
which are, in the mean, distributed uniformly in space and have time intervals with a Gaussian 
distribution. The generated clumps move radially outward with a velocity law determined by 
a f3 index, and the angular size of the clumps is assumed to keep fixed. By fitting the observed 
Cp/tTphot and the volume filling factor fv, the clump velocity law index /3 2) and clump 
, ejection rate Af 1) are inferred, and are found to be well constrained. In addition, the 

■ subpeak features on broad emission lines seem to support the clump ejection rate. Meanwhile, 

the fraction of the total mass loss rate that is contained in the clumps is obtained by fitting 
^ ■ the observed polarization. We conclude that this picture for the clump properties produces a 

I valuable diagnostic of WR wind structure. 
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H ■ 1 INTRODUCTION 



The intrinsic polarization of hot stars results from the scattering of starlight by electrons in aspherical 
stellar winds. There is evidence that the asymetry is not from a fixed structure but rather from stochastic 
ejection of clumps from the base of hot star winds. This evidence has been accumulating from a variety 
of spectropolarimetric studies that show that variability appears at a wide range of timescales of days, 
weeks, and months (Lupie & Nordsieck 1987; Taylor et al.l991). In the case of Wolf-Rayet stars, detailed 
observational data on the photometric, polarimetric, and spectral line profile variability have been presented 
and discussed by numerous authors (e.g., St. -Louis et al.l987; Drissen et al.l987; Robert et al.l989; Drissen 
et al.l992; Moffat & Robert 1992; Robert 1992; Lepine et al. 2000). Moffat & Robert (1992) and Lepine 
& Moffat (1999) have interpreted the Wolf-Rayet observations in terms of a distribution of dense clumps 
which result from hierarchical turbulence in the stellar winds. Brown et al. (1995) have discussed some of 
the physical properties of the larger clumps which dominate aspects of the data. Brown (1993) and Brown 
et al. (1999) conclude that clumps arise either by localized mass loss enhancements at the stellar surface or 
by the action of radiatively driven shocks sweeping up material on large scales. Brown et al. (2004) have 
discussed the combining effects of clumping and multiple scattering on the "momentum paradox" of WR 
star winds. 

With the systematic monitoring campaigns for WR stars, Robert et al. (1989) and Robert (1992) find 
that statistically, WR stars display polarization about 0.1%, and show broad band polarimetric variation 
(dp w 0.1% — 0.02%) that is much smaller than the fractional photometric variability (cTphot), the mean 
ratio being about TZ = (Tp/cTphot ~ 0.05. Richardson et al. (1996) investigate the statistical effect of having 
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larger numbers of clumps present and conclude that the clumps must be very dense so that their emission 
(cx n^V) is large enough to increase (iphot and/or their optical depth is large enough to reduce Cp by 
multiple scattering. Li et al. (2000) revisited the analysis by carrying out numerical simulations, following 
clump flow from the star with a f3 velocity law while accounting for the occultation of clumps behind 
the stellar disk, but while also retaining the single scattering assumption. They concluded that the clumps 
follow a large (3 velocity index, which means that the clumps are accelerating relatively slowly with radius 
in the wind. They also derived a certain range of clump ejection rate that along with /? allows for a fit to 
the statistical data. Valuable constraints on the WR parameters, such as mass loss rate, can be inferred from 
the observed mean polarization p (about 0.1% level), the polarimetric and photometric variances. These 
conclusions are supported when considered in conjunction with the number of distinct narrow features that 
are seen in emission line profiles. 

As for the observational evidence on clumping of hot stars in other wavebands, Ignace, Quigley, & 
Cassinelli (2003) used the observations with the Infrared Space Observatory (ISO) SWS spectrometer to 
constrain the velocity law and wind clumping, and found that a /3-value of the velocity law of the wind 
is about 2-3 and the volume filling factor is up to 40% for WR 136. Marchenko et al. (2006) obtain 
numerous overdense clumps in the wind of WR 135 with FUSE observations. Chandra high resolution 
spectral observations of line profiles have shown the clumping wind properties of hot stars (Owocki & 
Cohen 2006; Oskinova, Hamann, & Feldmeier 2007; Cassinelli et al. 2008). 

Recently, Davies et al. (2007) used the clump-ejection model to study the polarimetric variability of hot 
stars, in particular of luminous blue variables (LB Vs), and got the detailed clump parameters through their 
simulations. One of their conclusions is that many tiny clumps are ejected around hot stars. Of particular 
interest to us is that Davies et al. (2007) find a flaw in assumptions of Li et al. (2000), in their assumption that 
both the thickness of a clump is fixed and that the number of electrons in a clump is constant. These results 
from Davies et al. (2007) motivate us to update the model of Li et al. (2000). It is our goal to understand 
the general nature of the Wolf-Rayet stars, using a picture for the ejection of clumps into the wind. Here we 
allow for the thickness of a clump to vary with the wind flow. We find that the main conclusions in Li et al. 
(2000) still hold, and that the present model leads to well constraints on the parameters, such as [3 and Af. 

The basic model scenario is presented in Section 2, where the basic formulation concerned with the 
polarization and scattering light of both a single clump and ensemble of clumps are taken into account for 
various (3 and Af. In Section 3, the model results are presented and discussed, and in Section 4 general 
conclusions are presented. 



2 CLUMPING MODEL 

2.1 Basic Formulation of Polarization for A Single Clump 

We assume that the polarimetric and photometric variability of a WR star is due to localized mass ejections 
at the stellar surface, over which the clumps are generated at random positions, and at random time intervals 
with a normal distribution of mean value At. The clumps are then taken to move radially outward with a 
velocity law, thus the thickness Ar of the clump should simultaneously expand along the radial wind flows, 
while the solid angle Aft is assumed constant (The geometry of a single outflowing clump is shown in 
Fig.[T]). The electron density then decreases as r^^ v^^{r). The clumps thus have axisymmetric shapes and, 
on the assumption that they are not optically thick in the continuum, the results of Brown 8c Mclean (1977) 
can be used to find the polarization of a single clump as: 

P = Topt (1 - 37) sin^ i, (1) 

where 

r2 /-Ml 



Topt = y^ctt J J n{r, ^i) dr dfi (2) 



is a mean optical depth, and 



is a "shape" factor, i is the clump axis inclination to the line of sight; /i = cos i?, where is the clump 
opening angle between the axis of symmetry and the direction of the scatter seen from the center of the star; 
n{r, ji) is the electron number density in the clump; and (Tt is the Thomson scattering cross-section. For 
the local reference frame (r, -d) chosen, we let jix = \. Given constant solid angle AVl, the electron density 
is assumed to vary only with the distance, i.e., as n(r, ji) = n{r). To calculate the electron density in one 
clump, we use the mass conservation law: 
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where Mb is the mass outflow rate within one clump, p{r) is the mass density, and v{r) is the radial velocity 
law that the clump will follow, which we adopt to be of the common form: 

v{r)^v^ (5) 

where Voo is flie terminal clump speed, i?* is the photospheric radius of the WR star (We note that there are 
dynamical and effective optical photosphere notations, see Brown et al. 1995), 6 is a dimensionless param- 
eter to ensure that the initial wind speed is non-zero (6 — 0.995 is adopted throughout the simulations), and 
/3 is a velocity law index, one of our basic clump parameters. 
The electron density in a clump thus becomes 

P 

Ue — 

Me mH 

He niH Afi V 



fie rnn Ail v ' 

where niH is the hydrogen mass, /ig is the mean particle weight per free electron. 
To deal with the finite star geometry, the point source depolarization correction factors 
D{r/Rt), C{r/Rf,x) can be employed according to Cassinelli et al. (1987) and Brown et al. (1989), 



D 



-\ (7) 
and 

^ 8- J(l + i^)(l-3cos^x) .o. 

3(l + i?)(l + cos2x) ' ^' 

where x ~ r/R^,, again i?* is the photospheric radius of the WR star, and x is the scattering angle. 

We wish to combine Eqs.[T]to[8]to yield an expression for the polarization from a single clump with 
the assumed geometry. After the dimensionless treatment of r (i.e., rl/i?* = x\ and r2/i?* = x2), the 
polarization expression now becomes 



3 f X 

3 

= 77 '^T '^o R* (1 - M2)(M2 + M2) sin^ X 
Id 



X \ D{x)dx 



^2 / ^ 



xi 



a; \^ 1 / 1 



^x ~ b / x^ \ X 

and the scattered light intensity /s as a fraction of L^/Att in terms of Brown et al. (1995) is: 



fs = 



3 /"^' 



X \ C{r,x)dx 



X — b 



8- D{1 + D)il-3cos^ X) 



X — b 



dx. 



(10) 



4 



Q. Li, J.P. Cassinelli, J.C. Brown & R. Ignace 



2.2 Inference of Conservation of Electrons in A Clump 

It is assumed that the clump has an initial extent of radial thickness Ax — x2 — xl (= 0.01, for example) 
and solid angle Aft (= 0.04, for example), and the solid angle remains fixed. As the clump moves outward 
radially, obeying the velocity law, its thickness will naturally change to expand radially with both outer face 
and inner face obeying the local velocity expression. We may explicitly obtain the number of electrons in 
one clump 



p27r pfi2 nr2 

^ / / n,dV 

Jo Jfj.! Jrl 

^ MbR. dx ^^^^ 

where dV — —r^dfi d(f)dr in a local spherical coordinate, and for the geometry of the clump assumed, one 
may get = 1 and Ail = 2tt{1 — fi2). 

In Eq.[TTl for given values of f3 and initial Ax, the integral can be solved analytically and/or numerically. 
In the work, we adopt three acceleration cases: /3 — 0.5 for a rapidly accelerating flow; /3 = 1, a commonly 
assumed value for hot star winds; and (3 = 2 for the slow acceleration case. These are chosen so we can 
determine which case provides an improved fit to the WR observational properties. 

(1) In case of /3 = 0.5, we may analytically do the integral in the Eq.fTTIand obtain 



MbR* 



IJ-e rriH Voo 

MbR* 



x2 



1=1 Vx -b 



dx 



\/ x{x — b) + b In + ^/ x — b^ + conl 



IJ-e rUH Voo 

= Neo I ^2^2(3:2 - b) - y^xl{xl - b) + 61n 



x2 

xl 



x2 + \Jx2 - b 



= N, 



eO 



\/ x{x — b) + 6 In (^\/x + \/x — b^ + conl 



xl + ^/ xl 

1+0.01 



(12) 



where conl is a constant from integral and A^eo 



is also a constant for given the stellar param- 



/ie "Iff Voo 

eters. Owing to the boundary condition assumed, i.e. x2 = xl + Ax = 1 + 0.01 at .tI = 1, the expression 

1+0.01 

equals to 0. 103. Therefore, Eq.[T2]can clearly yield the 

1 



•y/ x{x — b) + b In {y/x + \/a; — 6) + conl 
following equation 



^x2{x2 -b)- y^T(^r~6)+61n^:^i^^S = 0.103. 

V xl + \' xl — b 

(2) In case of /3 = 1, we use the same method above and get the equation 



a:2 — xl + 6 In ■ 



x2-b 



1.103. 



xl-b 

(3) In case of /3 = 2, we repeat the same process and obtain the following equation 



x2- xl + 2b\n 



x2-b P{x2-xl) 
xl-b ^ {xl - b){x2 - b) 



= 134.2. 



(13) 



(14) 



(15) 



Eq.[T3Jll4l or[T3]sets up the coherent relationship between xl and x2. We may use them with a given 
value of (3 to get x2 once xl is specified. In the simulations, we use the Newtonian bisection and bracketing 
methods to solve the Eq. [T3] [14] or [T5] due to the nonlinear relationship between xl and x2. Note, it is 
clearly shown that the number of electrons in a clump with various (3 are different {N^ = A"eo x 0.103 for 
(3 = 0.5, Ne = Neo X 1.103 for /3 = 1 and Ne = A^eo x 134.2 for /3 = 2). [We may also have a special 
case of /? = 0, which will give the value of the integral equals to Ax (Aa;=0.01) and then we will have the 
expression A'e = A"eO x 0.01 for (3 = 0.] The electron scattering in the clump results in the polarization. 
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2.3 Relation between Time t since Clump Expulsion and the Clump Distance r from the Star 

The radial thickness varies with time as the clump moves out. We can connect the time since expulsion of a 
particular clump with its current radial distance together from the velocity law. Similarly, given time t, the 
clump distance (say r) can be obtained. From the clump velocity law, this relationship is. 



dr 

v{r) 



dx 



^ 7 (1 - b/xy 

dx 



[i^b/xY' 



(16) 



where t = i?*/uoo is denoted as the "flow time scale". 
Again, in the case of (3 — 0.5, 



^/x — b 



zdx 



^/x{x - b) - - In 



- 1 



con2. 



(17) 



If we adopt x = 1 sit = Q and use b — 0.995 so far, then we will get con2 = —0.14 in the Eq. [TtI 
Therefore, the expression between time and distance is 



- 1)2] - 0.14. 



In the cases of /3 = 1 and (3 — 2, the expression between time and distance are, respectively, 

t 



- = x + bln{x - b) + A.27 

T 



(18) 



(19) 



and 



t b 

- = x-b + 2b ln{x ~ b) 

r X — b 



208.55. 



(20) 



Therefore, from the numerical experiments, t/r may be stepwise accumulated by doing summation of 
At/r, and At/r is chosen to have a Gaussian distribution. 



2.4 Polarization for An Ensemble of Clumps 

The expressions above for polarization and scattering are for a single clump. For accounting for many 
clumps, we use the same approach as in Li et al. (2000). We set up a spherical coordinate system (r, 9, 4>) 
with the polar axis oz being along the line of sight. Then for each clump the "inclination" angle i is identical 
to both the scattering angle x ™d the polar angle 9, while the polarization position angle in the sky, ip, is 
just the coordinate component (j) (see Fig.[T]i. 

For a system of clumps labeled j = 1, A^, the total scattered light fraction /g and the net polarization 
are as usual for the optically thin case, simply given by the sum over j the Stokes intensity parameters 
Qj = pj cos20j, Uj = pj sm2(j)j of each to get the totals of Q and U, then finding p = {Q^ + U^y/"^ 
and position angle ^ = | arctan ^. It is to be understood that summations exclude all occulted clumps 
- i.e. clumps whose coordinates Xj, 6j, (f)j satisfy 9j > tt/2 and xj sm6j < 1. Note that fij = cos^j is 
uniformly sampled in the interval —1 to +1 and from to 27r . The radii Xj are determined by time and 
the velocity law. 

The total mass loss rate M is distributed among all clumps and the ambient, inhomogeneous but, on 
average, spherical "wind". We choose the fraction of the mass loss in the clumps to be 77. Hence, rjM = 
N N^li^rriH, where N is the mean clump ejection rate (s^^) and again is the total number of electrons 
in one clump. Thus, if we fix M and increase N then there are more clumps in any given range of r but each 
would be of smaller N^- We denote Af = N t, the number of clumps ejected per characteristic flow time 
(r = R^/voo), as the measure of the clump ejection rate. If we assume the portion of the total mass loss rate 
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7yM = kMb- Since in a flow time r, there are Af clumps ejected, we can attain k ~ Af. Therefore, we may 
replace Mi, in Eq.|6]with r/M/Af when accounting for many clumps in the simulations. 

The system of equations governing the time varying polarization and scattered light is then 



and 



N 



N 



ctt ^* (1 ~ M2) (M2 + M2) sin^ COS 24) j 



16 



3 n' ^ 



2:2 / ^, \ /3 



\- 1 1 



2 1/1 dxj 



.2/1 o^ / -"^ 



Y^s\v?e, / -3 /l - -jdx, (21) 



N 



3 = 1 

N 



3 ^ 

= t^ctt Uo R*{1 - M2)(/^2 + ^ sin^ 9j sin2(/)j 



16 



3 n' ^ 
= Ye^'^]^^*'-"^ ~ /^2)(M2 + ^4)z2 ^^""-^ sin20j 



^2 / ^. \/3 



X, 1 1 



2^/1 2 d'^j 

^ I-X2 / ^ \0 



= |gsin^.,.,i„2*£(-ii^) |^l-;|<i.,. (22) 

where nj, = -qM / AJ7 i?,^ Wqo) and 

3 



PO = 7^crTn-oi?*(l - M2)(M2 + M2) 

16 



3 r;i\f 
16 

Hence, the total polarization is given 



— brTTT— ^*(l-M2)(Ai2+M2)- (23) 



p = VQ^+TP (24) 

and 

N 

fs ^ ^. /si 



WR wind clump polarization variability 



7 



where 



3 

—aTnoR*{l - /U2) 



X 



AA^A, 3(1 + " ^ ^ 



3 

/o = — crT?^oi?*(l - M2) 



^777^ R,{l-^l2). (26) 



3 r;Af 

— T -I 

16 /leniH An R,: Vo 

Note again, in the above expressions, both xl and x2 are varying with time, so we need to use Eq.fTSi 
[l4l or[T5]in regards to various /3 to determine x2 once xl is given. Here, xl is obtained from Eq.[T8][T9] or 
l20lin case of different f3, if time t is known from the numerical experiments. 

On inspection of the above equations for polarization and scattered light, we expect that, for given M, 
AO, initial A?-, and /3, results should depend mainly on J\f. For a fixed mass loss rate and flow time, if the 
clump generation rate is low, only a few clumps each of large density will be present near the star and these 
dominate the p and fs values, as is shown in Li et al. (2000). But for high generation rates, many low density 
clumps near the star will be controlling p and fs- So the same total number of electrons is redistributed in 
different number of clumps, resulting in different statistical means and variances in the polarization and 
scattered light fraction. For a fixed clump ejection and mass loss rate, the number of clumps in the inner 
radii near the star has a steady mean value and so therefore do the resulting polarization, scattered intensity 
and their variances, but these values change with Af, AI, and /3. So their observed values allow inference of 
the clump emission and flow parameters. 



3 MODEL RESULTS AND DISCUSSION 

1. As a start, let us calculate the thickness, volume and polarization of a single clump with varying clump 
location, supposing sinx = 1 and /3 ~ 1. The results are displayed in Fig|2] In Fig. 2a the solid line 
denotes the inner radial boundary (say xl) and the dotted line denotes the outer radial boundary (say 
x2). We see that as the time progresses, x2 increases faster than xl, which directly induces the thickness 
expansion. In Fig.2a, as for comparison, x2 {x2 = xl + Ax) with the constant thickness Ax = 0.01 is 
shown in the dashed line which almost overlaps the solid line xl. In Fig. 2b, the solid line denotes the 
case for which a constant clump thickness is assumed. In contrast, the dotted line shows the expansion 
of the thickness with its location (i.e. x), and we see that the thickness increases rapidly among a; = 1 to 
X = 2, where the winds are mainly accelerated. In Fig. 2c, the solid line denotes the volume of the clump 
with constant thickness Ax = 0.01 and the dotted line denotes the volume of the clump with expanding 
thickness, which is increasing dramatically faster than the constant thickness case. We compare the 
polarization of varying thickness (the dotted line in Fig. 2d) with that in the Davies et al. (the solid line 
in Fig. 2d) for one single clump and see that in both cases, the polarization decreases as the clump moves 
outward, while there is a peak polarization about a; = 1.2 for the latter case. 

2. Furthermore, we calculate the thickness, volume and polarization of a single clump with varying clump 
location, supposing sinx = 1 but with different jS, as is shown in Fig. 3. To see how the xl varies with 
time (say i/r) in case of differing /3s, we show the results in Fig. 3a in which fast acceleration of /3 = 0.5 
is displayed in dashed line and moderate acceleration of /3 = 1 is in solid line, and slow acceleration 
of = 2 is in dash-dotted Une. We also show the radial thickness varies with location of the clump in 
Fig. 3b. The dash-dotted line denotes the case of /3 = 2, the solid line denotes the case of /3 = 1, and the 
dashed line denotes the case of f] — 0.5. In addition, we show the constant thickness case with (3=1 
in dotted line. In Fig. 3c, we show that the volume changes with the clump location in various cases of 
f3 as denoted in the figure. In Fig. 3d, we show the polarization varies with location of the clump. The 
dash-dotted line denotes the case of /3 = 2, the soUd Une denotes the case of /3 = 1, and the dashed Une 
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the start values of polarization in various case of (3 are different since the electron number in one single 
clump for various case of (3 is different, as is shown from Eq.[T2]which is related to /3. And polarization 
is mainly determined by the number of electrons in the clump (Brown et al. 1995). 
3. Taking into account the behavior of an ensemble of clumps, we compute the polarization, the fraction 
of scattering light intensity, and their variation. In particular, we calculate the ratio of their variations 
for various clump ejection rate A/" in a flow time with various (3 = 0.5, 1, 2. Interestingly, <Tp/p and 
Cp/cphot do not depend on any specific star, but only on the JV and (3. It is not surprising at this since 
p (or /s) linearly relies on po (or /o) which is determined by the stellar parameters of a star So we 
enable to treat (Tp/cTphot as a probe to explore the probable Af and (3. To get rid of the influence of any 
specific star, we divide p (as well as Q and U) in Eq.|24]by po in Eq.|23]and divide fs in Eq.|25]by /o 
in Eq.|26] The results are displayed in Tables[Tl|2]and[3] Using the values in these tables, we plot (Jp/p 
and CTp/tTphot with various J\f in cases of various /3 — 0.5, 1, 2. The observed value of cTp/p — 0.5 
is denoted in dotted line in the upper panel of Fig. |4] and the observed value of (Tp/o-phot ~ 0.05 is 
denoted in dotted line in the lower panel of Fig.lD We found that to achieve the observed value, three 
patterns could hold, for /3 = 0.5 — » TV - 500; for /3 = 1 — > TV - 150; and for /3 = 2 — > TV - 1. 
Although there seems no solid unique solution for (3 and TV in terms of the above considerations, while 
combining the effects of TV and /3 on CTp/uphot and fv the volume filling factor, it seems the observed 
data might prefer the case of TV = 1 and (3 — 2. We will further demonstrate this later 
Hamann and Koesterke (1998) found the clump volume filling factor of WR subtype WN stars about 
30% in terms of their spectrum analyses. In the work of the paper, we are able to calculate the clump 
volume filling factor from our model. After doing summations of volumes that all clumps occupy and 
divided it by the whole space accounted for, then the volume filling factor denoted as fy, can in common 
definition be obtained. 

We present the model results of fv in column 9 of Tables [T] |2] and [5] in case of various /3s. We also 
plot fv vs TV in Fig. |5] in which the dotted line denotes the "observed" of 30%. The figure shows 
/3 = 0.5 — > TV - 700; /3 = 1 — > TV - 100; (3^2 — > TV - 1, in regards to the observed. Note, 
fv greater than 1 means clumps overlapped or merged. So for the assumption of thin clumps, we would 
rule out the cases of fv > 1- 

In general, the hot star winds show two-component property: the smooth ambient and clumps. Of course, 
when the clumps move outward to very far away distance, they will become as part of the interstellar 
medium. However, the density of clumps is in fact higher than that of the ambient wind around the star 
within several hundred of stellar radius. The smaller the value of the volume filling factor, the stronger 
the wind clumping. A clumping wind may cause higher emission for the same amount of material, which 
implies that the mass loss rates by spectroscopic analysis under the assumption of a smooth wind are 
systematically overestimated by typically a factor of ~ 2 and even higher, and the mass loss rates 

V fv 



in turn consequently affect the stellar structure and evolution. Such clumping presumably arises from the 
inherent instability of radiation driven winds and would influence the strength of the electron scattering 
wings, change ionization and line ratios, and cause polarization variability and profile variability. In 
fact, Abbott et al. (1981) has taken the inhomogeneities (i.e. clumps) in the wind into account and got 
the relationship of the mass loss rate, radio flux, and the filling factor Later, in most of models, in 
most simplicity, the inter-clump medium is usually treated as void, which recently leads to the porous 
scenario (Owocki and Cohen, 2006) for OB stars, as the photons are able to leak freely through the large 
separation between clumps. However, our model clearly indicates that the void treatment for WR stars 
might likely cast some doubts, since the inter-clump medium, occupying a very large mass percentage 
of the total, may play a role in the hot wind emission, in particular, near the star. The whole inter-clump 
medium, being of the smooth spherical wind ambient, globally contributes nil to the net polarization 
due to the cancellation effect of polarization. 

We compute the photometric and polarimetric intensity and their fluctuations in case of various (3s. 
Figs.l6l[8]and[T0lshow how polarization, position angle, and scattered light intensity change with the 
total number of clumps emitted from the start. Polarization changes are also shown as a locus in the Q - 
U plane (see Figs.|6l[8]and[T0li which, as expected, shows no strong preferred direction, since the mean 
structure is quasi-spherical. In Figs.|7]|9l and[TT]we show "observational" time-smoothed results for the 
variations in mean polarization and scattered light. Standard deviations, a, of these quantities are also 
plotted in these figures. The ratio TZ versus total number of clumps TV as time progresses is plotted in 
Fig. [El from which we see that in order to obtain the sustainable TZ = 0.05 as observed, the pattern of 
/3 = 2 and TV = 1 is indeed preferred. This seems comparable to the number of subpeaks of the broad 
Une observed. Given the model quantities i?* = IO-R0 and Wqo — 1800 km s^^, the flow time r is about 
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new clumps will occur. Hence, there would consequently result in 10 more detectable subpeaks since the 
clump emission is dominated by the inner clumps in the line emission regions (LER), which are closer 
to the star This conclusion is indeed consistent with the observations and their analysis (Robert 1992, 
Brown et al. 1995). Note, envisaging Tables[Tl|2]and[3]in which TZ = 0.05 seem to be attained, however, 
observing Figs.|6l[8][T0l and their time averaged results in Figs.|2l|9][TT] in particular[l2l we found that 
only the case of /3 = 2 and J\f = 1 shows a real and sustainable result. Neither the case of /3 = 1 and 
Af = 200, nor the case /3 — 0.5 and Af — 1000, could achieve a long-standing TZ = 0.05. Therefore, 
we have to abandon these cases. If we plot the number of clumps versus the distance (Fig.fTSI). we will 
find that in case of various (3 — 0.5, 1, 2 but with same Af — 1, there will several hundred of clumps 
radially staggering outward from li?* to 5i?* in case of jS — 2 being very slowly accelerated, but just a 
few in case of = 0.5, 1 being rapidly accelerated. In terms of the concept of LER in Lepine & Moffat 
(1999) and Dessart & Owocki (2005), the LER with velocity space is about among 0.4-0. 9tioo where 
the clumps are accelerated. We note that the WR-wind acceleration length scale /3i?* ^ 20Rq in their 
results, is compatible with our results, with i?,* = lOi?© applied and /3 = 2 inferred. 
5. For the fraction factor of the mass loss rate into the clumps 77 in regards to po in Eq. |23]and /o in 
Eq. |26] we could use the typical polarization observed for WR stars p = 0.1 — 1% to constrain it. 
If the following typical parameters are employed: M — 2.5 x lO~^M0/year, Afi = 0.04, initial 
Aro = O.Oli?*, fie = 2, i?* = lOi?0 and Woo = 1800kms"\ then we may gain po = 0.01477^. 
From the above discussion, the pattern (3 = 2 and Af = 1 is preferred, then we check Table [3] and get 

p' = p/po ^ 124.85 — > p = 1.8477. Therefore, we enable to gain 77 = p/1.84 ~ lO^'' 10^^. This 

would imply that only a quite small portion of the winds is going to deposit into the clumps but most 
of winds are blown as an ambient where clumps embedded. In contrary, one might recall the solar wind 
picture in which the solar wind is a coherent outward expansion of the solar corona, frequently involving 
the corona mass ejection (CME) events. Hence, it is not surprising that there are two components in the 
hot star environments, being the ambient winds and the clumps (or wind-blown bubbles). However, 
the formation mechanisms of the solar winds and hot star winds are distinctly different. The former 
is driven by the gas pressure gradient of the high temperature solar corona, and the latter is driven by 
the pressure of the radiation emitted by the hot star, so-called the continuum-driven and line-driven. 
Note that the physical quantities in Tables [T] |2] and |3] such as p' = p/po, a'^ = o-p/po, /g = /s//o, 

"'phot = o-phot / fo, vl fs, CTp/fTphot, scale with pQ or /o, so they are dimensionless numbers. The fraction 
factor 77 only affects the values of po and /q. Thus for a specific star one could gain 77 incorporating 
the stellar parameters and observed polarization. We realize that the clumps ejected in our model are 
massive- or macro-clumps, distinguished from the local perturbations in the winds. 

4 CONCLUSIONS 

In this paper, we update the previous model proposed by Li et al. (2000) and apply the clump ejection sce- 
nario to explain the WR star wind polarization and its variability observed, by accounting for the expansion 
of the clumps along the wind flows but keeping the solid angle constant. We may gain the main conclusions 
as follows. 

(1) From numerous model simulations using various /3 and Af, incorporating with the volume filling 
factor fv, we found /? ~ 2 and J\f ^ 1 are preferred for explaining the observational data on photometric, 
polarimetric, and spectral line profile variability of WR stars. This gives well model constraints on the stellar 
wind properties of hot stars. 

(2) We also found that a small fraction as 10^^ of the wind material deposits into the clumps but most of 
the mass loss are going to the space as wind ambient. This quantitative estimation of the mass fraction into 
the wind ambient implies that the porous wind models for WR stars might be cautious since the inter-clump 
medium is far from void. 

In summary, this updated model which is improved from Li et al. 2000, with inclusion of the radial 
expansion of the thickness of clumps but retaining the dependence on the /? velocity law and stellar oc- 
cultation effects, yields results consistent with observations in regards to the mean polarization p, the ratio 
TZ = CTp/cTphot of polarimetric to photometric variability, and the volume filling factor fv - It offers a quanti- 
tative estimation of the mass fraction into the clumps and ambient winds as well. Hence, the model produces 
a valuable diagnostic of WR wind structure. 
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Table 1 Simulation results for finite star source with occupation and velocity law with /3 = 0.5 



N 


u 


P' 




/a 


"■phot 


p/fs 


fp/cTphot 


fv 


5000 


1/4 


0.8394 


0.8632 


4.928 


4.470 


0.337 


0.382 


0.0001 


5000 


1 


0.3565 


0.2032 


2.319 


1.138 


0.305 


0.353 


0.0003 


5000 


5 


0.1557 


0.0816 


1.927 


0.276 


0.160 


0.479 


0.0016 


5000 


10 


0.1091 


0.0557 


1.891 


0.234 


0.114 


0.473 


0.0033 


5000 


50 


0.0500 


0.0265 


1.835 


0.195 


0.054 


0.271 


0.0164 


5000 


100 


0.0364 


0.0204 


1.783 


0.239 


0.040 


0.168 


0.0328 


5000 


200 


0.0268 


0.0146 


1.693 


0.304 


0.031 


0.094 


0.0654 


5000 


1000 


0.0143 


0.0051 


1.229 


0.424 


0.023 


0.024 


0.3190 


5000 


2000 


0.0102 


0.0031 


0.910 


0.392 


0.022 


0.016 


0.6239 



Notes: Column 1 is the total number of clumps applied in the simulations. Column 2 is the clump ejection rate in a flow 
time. Colurrm 3 is p' = p/po. Column 4 is Up = a-p/po. Column 5 is = /s//o. Column 6 is CTphot = (Jphot/ fo- 
Column 7 is p/fs. Column 8 is Opjaphot- Column 9 is the volume filling factor fv- 
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Table 2 Simulation results for finite star source with occultation and velocity law with (3—1 
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N 


JV 


P' 


/ 

CTp 




/ 

"^phot 


P/fs 


fp/Cphot 


fv 


5000 


1/4 


i.mi 


2.1003 


22.762 


17:995 


0.252 


0.231 


0.0015 


5000 


1 


1.4512 


0.7755 


15.340 


5.933 


0.187 


0.259 


0.0035 


5000 


5 


0.6270 


0.3205 


14.144 


2.371 


0.088 


0.268 


0.0175 


5000 


10 


0.4432 


0.2223 


14.077 


1.753 


0.062 


0.251 


0.0350 


5000 


50 


0.2066 


0.1172 


13.805 


1.664 


0.030 


0.140 


0.1725 


5000 


80 


0.1641 


0.0982 


13.565 


1.985 


0.024 


0.098 


0.2727 


5000 


100 


0.1458 


0.0904 


13.395 


2.170 


0.022 


0.082 


0.3381 


5000 


150 


0.1189 


0.0703 


13.021 


2.557 


0.018 


0.054 


0.4962 


5000 


200 


0.1054 


0.0510 


12.663 


2.852 


0.016 


0.035 


0.6460 


5000 


400 


0.0788 


0.0279 


11.316 


3.529 


0.014 


0.016 


1.1568 


5000 


1000 


0.0506 


0.0153 


8.065 


3.647 


0.013 


0.008 


2.2168 


5000 


2000 


0.0255 


0.0088 


4.986 


2.591 


0.010 


0.007 


4.0601 



Notes: Parameters have the same meaning as that in Table [T] 



Table 3 Simulation results for finite star source with occultation and velocity law with (3 = 2 



N 


N 


-/ 
P 


/ 

CTp 




/ 

''"phot 


P//b 


C"p/o"phot 


fv 


5000 


1/4 


247.339 


126.392 


17517.453 


3157.450 


0.028 


0.079 


0.1071 


5000 


1 


124.8503 


64.1338 


17300.174 


1918.061 


0.014 


0.066 


0.4098 


5000 


5 


48.1358 


28.5260 


16625.262 


2933.660 


0.006 


0.019 


1.7094 


5000 


10 


33.3430 


17.1025 


15757.138 


3900.831 


0.004 


0.009 


2.6139 


5000 


50 


14.3191 


3.7360 


9381.283 


4799.212 


0.003 


0.002 


1.3993 


5000 


80 


9.6549 


2.6050 


6372.940 


3492.983 


0.003 


0.002 


3.3746 


5000 


100 


7.8479 


2.1951 


5201.766 


2891.146 


0.003 


0.002 


5.6845 


5000 


150 


5.2524 


1.4971 


3550.873 


2009.073 


0.003 


0.002 


12.349 


5000 


200 


3.8611 


1.0675 


2689.970 


1536.053 


0.003 


0.001 


11.362 


5000 


400 


1.9001 


0.5064 


1352.113 


778.007 


0.003 


0.001 


15.412 


5000 


1000 


0.7599 


0.2026 


540.961 


311.203 


0.008 


0.001 


15.406 


5000 


2000 


0.3800 


0.1013 


270.565 


155.602 


0.003 


0.001 


15.396 



Notes: Parameters have the same meaning as that in Table[T] 
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Fig. 1 Geometry of scattering from one clump. We employ a spherical coordinate system 
{r, 9, (p) with oz along the line of sight. For each clump the "inclination" angle i is identical 
to the scattering angle x as well as the polar angle 9, while the polarization position angle on the 
sky, ip, is the coordinate component (p (see Li et al. 2000). 
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Fig. 2 Model results are shown for a single clump: (a) The radial extent vs time. The sohd line 
denotes the inner radial boundary (say xl) and the dotted line denotes the outer radial boundary, 

x2 (x2 = xl + Ax) with the constant thickness Ax = 0.01, is shown in the dashed line which is 
almost overlap the solid line xl. (b) The clump thickness versus the clump location. The solid line 
denotes the constant clump thickness as assumed. In contrast, the dotted line shows the expansion 
of the thickness with its location (i.e. x). (c) The clump volume vs its location. The sohd line 
denotes the volume of the clump with constant thickness Ax = 0.01 and the dotted line denotes 
the volume of the clump with expanding thickness, (d) Polarization vs the clump location. The 
dotted hue is for the varying thickness and the sohd hue is for the constant thickness as in Davies 
et al (2007). 




Fig. 3 Model results are displayed in various cases of /3 as denoted in the figure for one single 
clump. 




Fig. 4 The upper panel is ap/p vs Af in various cases of /? as denoted in the plot. The lower 
panel is Cp/crphot vs Af. The observed are denoted in dotted lines in two panels. 
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Fig. 5 Volume filling factor /v vs M in various cases of /3 as denoted in the plot. 
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Fig. 6 The figures show respectively, instant by instant, model results vs number N of clumps 
emitted thus far for the following: upper left panel - polarization p' = p/po; upper right - 
scattered light = /s//o; lower left panel - polarization position angle tp; lower right panel - 
Q' = Q{N)/po versus U' = U{N)/po. (3 = 0.5 and M = 1000 are applied. 
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Fig. 7 Based on the same data as Fig. |6] here we show smoothed results versus number of 
clumps N (increasing with time) for the following observables, with parameters as in Fig.|6l (a) 
mean polarization p' = p/po; (b) variance of polarization (t^ = <Jp/po', (c) mean scattered light 

fraction fs = /s//o; (d) variance of scattered hght a'p^^^^ = cTphot/fo- 



WR wind clump polarization variability 



19 



0.3 



0.2 



0.1 



I ' ' ' ' I 



_L i 1 I •) 



J L±J L 



5000 10^ 1.5x10* 2x10-^ 
N 




5000 10" 1.5x10" 2x10" 
N 



150 



100 




0.2 



b 



-0.2 



5000 



10^ 1 5^10" 2x10" 
N 




Fig. 8 The figures are similar to Fig. |6l except (3=1 and Af = 200 are applied. 
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Fig. 9 The figures are similar to Fig. |7] except /? = 1 and Af — 200 are applied. 
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Fig. 10 The figures are similar to Fig.|6l except [3 — 2 and J\f — 1 aie applied. 
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Fig. 11 The resulting figures are similar to Fig.|7] except (3 — 2 andAf = 1 are applied. 
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Fig. 12 Ratio {TZ = (jp/ o-phot) of polarimetric to photometric standard deviations vs number of 
clumps N. The upper panel (a) is for /3 = 0.5 and M = 1000. The middle panel (b) is for /3 = 1 
and Af = 200. The lower panel (c) is for /3 = 2 and M = I. The steady mean value TZ about 
0.05 can be sustainable in case of /3 = 2 and N = 1. 
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Fig. 13 The number of clumps N vs the radial distance (x) in case ofN= 1. The solid line is 
for /3 = 2, the dotted line is for /3 = 1, and the dashed line is for /3 = 0.5. Due to the "dwelUng" 
time around the star is longer in case of larger /3, the number of clumps are more, in comparison 
with the smaller (3 case. 



